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Abstract 

The prediction and prevention of traumatic brain injury is a very important aspect of preven- 
tive medical science. This paper proposes a new coupled loading-rate hypothesis for the traumatic 
brain injury (TBI), which states that the main cause of the TBI is an external Euclidean jolt, 
or SE(3)— jolt, an impulsive loading that strikes the head in several coupled degrees-of- freedom 
simultaneously. To show this, based on the previously defined covariant force law, we formulate 
the coupled Newton-Euler dynamics of brain's micro-motions within the cerebrospinal fluid and 
derive from it the coupled SE(3) — jolt dynamics. The SE(3) — jolt is a cause of the TBI in two 
forms of brain's rapid discontinuous deformations: translational dislocations and rotational discli- 
nations. Brain's dislocations and disclinations, caused by the SE(3)— jolt, are described using the 
Cosserat multipolar viscoelastic continuum brain model. 

Keywords: Traumatic brain injuries, coupled loading-rate hypothesis, Euclidean jolt, coupled 
Newton-Euler dynamics, brain's dislocations and disclinations 



1 Introduction 

Traumatic brain injury (TBI) continues to be a major health problem, with over 500,000 cases per 
year with a societal cost of approximately $85 billion in the US. Motor vehicle accidents are the leading 
cause of such injuries. In many cases of TBI widespread disruption of the axons occurs through a 
process known as diffuse axonal injury (DAI) or traumatic axonal injury (TAI) |Singh et al 2006 . TBI 
occurs when physical trauma causes brain damage, which can result from a closed head injur) 3 or a 
penetrating head injury^ In both cases, TBI is caused by rapid deformation of the brain, resulting in a 
cascade of pathological events and ultimately neuro-degeneration. Understanding how the biomechan- 
ics of brain deformation leads to tissue damage remains a considerable challenge [Morrison et al 2006] . 

Parts of the brain that can be damaged include the cerebral hemispheres, cerebellum, and brain 
stem. TBI can cause a host of physical, cognitive, emotional, and social effects |NIH 200 2, Rapp 2008 . 
Half of all TBIs are due to transportation accidents involving automobiles, motorcycles, bicycles, and 
pedestrians. These accidents are the major cause of TBI in people under age 75. For those aged 75 
and older, falls cause the majority of TBIs. Approximately 20 % of TBIs are due to violence, such as 
firearm assaults and child abuse, and about 3% are due to sports injuries. Fully half of TBI incidents 
involve alcohol use NIH 2002J . TBI is a frequent cause of major long-term disability in individuals 
surviving head injuries sustained in war zones. This is becoming an issue of growing concern in 
modern warfare in which rapid deployment of acute interventions are effective in saving the lives of 
combatants with significant head injuries. Traumatic brain injury has been identified as the 'signature 
injury' among wounded soldiers of the current military engagement in Iraq [Mason 2007, Hoge 2008] . 
Rapid deformation of brain matter caused by skull acceleration is most likely the cause of concussion, 
as well as more severe TBI. The inability to measure deformation directly has led to disagreement 
and confusion about the biomechanics of concussion and TBI Bayl y et al 2005| . 

TBI can be mild, moderate, or severe, depending on the extent of the damage to the brain. Out- 
come can be anything from complete recovery to permanent disability or death (see [Chen et al 2 008]). 



1 A closed injury occurs when the head suddenly and violently hits an object but the object does not break through 
the skull. 

2 A penetrating injury occurs when an object pierces the skull and enters brain tissue. 



1 



Some symptoms are evident immediately, while others do not surface until several days or weeks after 
the injury NI H 2002] , With mild TBI, the patient may remain conscious or may lose consciousness 
for a few seconds or minutes; the person may also feel dazed or not like him- or herself for several 
days or weeks after the initial injury; other symptoms include: headache, mental confusion, light- 
headedness, dizziness, double vision, blurred vision (or tired eyes), ringing in the ears, bad taste in 
the mouth, fatigue or lethargy, a change in sleep patterns, behavioral or mood changes, trouble with 
memory /concentration/calculation. With moderate or severe TBI, the patient may show these same 
symptoms, but may also have: loss of consciousness, personality change, a severe/persistent /worsening 
headache, repeated vomiting/nausea, seizures, inability to awaken, dilation (widening) of one or both 
pupils, slurred speech, weakness/numbness in the extremities, loss of coordination, increased con- 
fusion, restlessness/agitation; vomiting and neurological deficit together are important indicators of 
prognosis and their presence may warrant early CT scanning and neurosurgical intervention. 

In particular, standard medical statistics suggest that the loss of consciousness in boxing knock- 
outs and road-vehicle crashes is caused by rotation of the brain-stem Misra and Chakravart y 1984] 
as a dynamic response [Ivancevic and Ivancevic 2006a[ llvancevic and Ivancevic 2007dj of a head-neck 
system to an impulsive load |Misra and Chakravarty 1985| . It is generally associated to the fol- 
lowing three syndromes: Locked-In, Semi- Coma, and Akinetic Mute, all three characterized by the 
total loss of gesture, speech and movement. The cognitive abilities [Ivancevic and Ivancevic 2007c] 
Ivancevic and Aidman 2007 can still be intact, but the patient cannot express himself by either speech 
or gesture. Recall that the brain stem, including Midbrain, Pons and Medulla Oblongata, is located 
at the base of the brain. It is the connection between the cortex and the spinal cord, containing 
motor neural pathways for voluntary movement from the upper part of the brain. The brain stem 
also controls such automatic functions as breathing, heart rate, blood pressure, swallowing, sleep pat- 
terns and body temperature. Weaker injuries include another three symptoms: abnormal respiration 
(hyperventilation and abnormal breathing patterns: ataxic, clustered, hiccups); pupils: dilated, fixed; 
and movement (if any): abnormal extensor. 

The natural cushion that protects the brain from trauma is the cerebrospinal fluid (CSF). It 
resides within cranial and spinal cavities and moves in a pulsatile fashion to and from the cranial 
cavity (see FigureQ}. This motion can be measured by functional magnetic resonance imaging (fMRI, 
see [Spkoloff 2008] for a review) and may be of clinical importance in the diagnosis of several brain 
and spinal cord disorders such as hydrocephalus, Chiari malformation, and syringomyelia. It was 
found in [Maie r et al 1994] that brain and CSF of healthy volunteers exhibited periodic motion in 
the frequency range of normal heart rate. Both brain hemispheres showed periodic squeezing of the 
ventricles, with peak velocities up to 1 mm/sec followed by a slower recoil. Superimposed on the 
regular displacement of the brain stem was a slow, respiratory-related periodic shift of the neutral 
position. During the Valsalva maneuver, the brain stem showed initial caudal and subsequent cranial 
displacement of 2-3 mm. Coughing produced a short swing of CSF in the cephalic direction. The 
pressure gradient waveform of a linearized Navier-Stokes model of the pulsatile CSF flow was found in 
|Loth et al 2001] to be almost exclusively dependent on the flow waveform and cross-sectional area. 

The state of head injury biomechanics:: past, present, and future was presented in [Goldsmith 2001] , 
dealing with components and geometry of the human head, classification of head injuries, tolerance 
considerations, head motion and load characterization, experimental dynamic loading of human living 
and cadaver heads, dynamic loading of surrogate heads, and head injury mechanics. The subse- 
quent paper [Goldsmith and Monson 2005] described physical head injury experimentation involving 
animals (primarily primates), human cadavers, volunteers, and inanimate physical models. 

Besides motor accidents, concussion (or mild TBI), occurs in many activities, mostly as a result 
of the head being accelerated. For example, although a popular endeavor, boxing has fallen under 
increased scrutiny because of its association with TBI. The injury rate in professional boxing matches 
is high, particularly among male boxers. Superficial facial lacerations are the most common injury 
reported. Male boxers have a higher rate of knockout and technical knockouts than female box- 
ers [Bledsoe et al 2005] , Although the epidemiology and mechanics of concussion in sports have been 
investigated for many years, the biomechanical factors that contribute to mild TBI remain unclear be- 
cause of the difficulties in measuring impact events in the field. The objective of [Beckwith et al 2007] 
was to validate an instrumented boxing headgear (IBH) that can be used to measure impact sever- 
ity and location during play. Based on this study, the IBH is a valid system for measuring head 
acceleration and impact location that can be integrated into training and competition. Similarly, a 
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Figure 1: Human brain and its SE(3)-group of microscopic three-dimensional (3D) motions within 
the cerebrospinal fluid inside the cranial cavity. 



comprehensive study has been conducted by [Newman et al 2005] to understand better the mechanics 
of the impacts associated with concussion in American football, involving a sequence of techniques to 
analyze and reconstruct many different head impact scenarios. 

Impact biomechanics from boxing punches causing translational and rotational head acceleration 
was experimentally studied in [Viano et al 2005] . Olympic boxers threw four different punches at an 
instrumented Hybrid III dummy and responses were compared with laboratory-reconstructed NFL 
concussions. Instrumentation included translational and rotational head acceleration and neck loads 
in the dummy. Biaxial acceleration was measured in the boxer's hand to determine punch force. 
Hybrid III dummy head responses and FE brain modelling were compared to similarly determined 
responses from reconstructed concussions in professional NFL football players. The hook produced 
the highest change in hand velocity (11.0 +/- 3.4 m/s) and greatest punch force (4405 +/- 2318 N) 
with average neck load of 855 +/- 537 N. It caused head translational and rotational accelerations of 
71.2 +/- 32.2 g and 9306 +/- 4485 r/s 2 . These levels are consistent with those causing concussion in 
NFL impacts. However, the head injury criterion (HIC) for boxing punches was lower than for NFL 
concussions because of shorter duration acceleration. Boxers deliver punches with proportionately 
more rotational than translational acceleration than in football concussion. Boxing punches have a 
65 mm effective radius from the head eg, which is almost double the 34 mm in football. A smaller 
radius in football prevents the helmets from sliding off each other in a tackle. Similarly, impacts 
causing concussion in professional football were simulated in laboratory tests to determine collision 
mechanics |Viano et al 2 007 . This study focused on the biomechanics of concussion in the struck 
player, addressing head responses causing concussion in the NFL players. 

Head injury mechanisms are difficult to study experimentally due to the variety of impact condi- 
tions involved, as well as ethical issues, such as the use of human cadavers and animals 



Krabbel and Appel 1995 . A number of finite element (FE) models and analysis studies have been 
conducted in order to understand the mechanism of TBI. An FE analysis was carried out in [Chu et al 1994| 
to study the mechanism of cerebral contusion. Clinical findings indicate that most cerebral contusions 
in the absence of skull fracture occur at the frontal and temporal lobes. To explain these observations, 
cavitation and shear strain theories have long been advocated. Plane strain finite element models of 
a para-sagittal section of the human head were developed in the present study. The model was first 
validated against a set of experimental results from the literature. Frontal and occipital impacts 
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were then simulated, and pressure and shear stress distributions in the brain were compared. While 
comparable negative pressures always developed in the contrecoup regions, shear stress distributions 
remained nearly identical regardless of the impact direction, consistent with the clinically observed 
pattern for contusion. Therefore, shear strain theory appeared to account better for the clinical 
findings in cerebral contusion. 

A 3D FE model based on the anatomical features of the adult human cranium was developed in 
|Krabbel"an d Appc lT995| . The complex cranial geometry was measured from a series of 2D computer 
tomography images. The CT scans were transformed with a self-developed preprocessor into a finite 
element mesh. A review of the existing FE models was presented in [Voo et aT~ 1996 for the biomechan- 
ics of human head injury. More recent models incorporated anatomic details with higher precision. 
The cervical vertebral column and spinal cord were included. Model results had been more qualitative 
than quantitative owing to the lack of adequate experimental validation. Advances included transient 
stress distribution in the brain tissue, frequency responses, effects of boundary conditions, pressure 
release mechanism of the foramen magnum and the spinal cord, verification of rotation and cavitation 
theories of brain injury, and protective effects of helmets. These theoretical results provided a basic 
understanding of the internal biomechanical responses of the head under various dynamic loading 
conditions. The mechanism of brain contusion has been investigated in in |Huang et aH OOO using a 
series of 3D FE analyses. A head injury model was used to simulate forward and backward rotation 
around the upper cervical vertebra. Intracranial pressure and shear stress responses were calculated 
and compared. The results obtained with this model support the predictions of cavitation theory 
that a pressure gradient develops in the brain during indirect impact. Contrecoup pressure-time his- 
tories in the para-sagittal plane demonstrated that an indirect impact induced a smaller intracranial 
pressure (-53.7 kPa for backward rotation, and -65.5 kPa for forward rotation) than that caused by a 
direct impact. Comparison of brain responses between frontal and lateral impacts was performed in 
|Zhang et al 2001] by FE modelling. Identical impact and boundary conditions were used for both the 
frontal and lateral impact simulations. Intracranial pressure and localized shear stress distributions 
predicted from these impacts were analyzed. The model predicted higher positive pressures accom- 
panied by a relatively large localized skull deformation at the impact site from a lateral impact when 
compared to a frontal impact. Lateral impact also induced higher localized shear stress in the core 
regions of the brain. 

A nonlinear viscoelastic FE model for brain tissue was developed in Brands et al 2004] . To obtain 
sufficient numerical accuracy for modelling the nearly incompressible brain tissue, deviatoric and vol- 
umetric stress contributions were separated. Deviatoric stress was modelled in a nonlinear viscoelastic 
differential form. An attempt was made in | Zhang et al 2004] to delineate actual injury causation and 
establish a meaningful injury criterion through the use of the actual field accident data. Twenty-four 
head-to-head field collisions that occurred in professional football games were duplicated using a vali- 
dated FE human head model. The injury predictors and injury levels were analyzed based on resulting 
brain tissue responses and were correlated with the site and occurrence of mild TBI. Predictions indi- 
cated that the shear stress around the brainstem region could be an injury predictor for concussion. 
The controlled cortical impact model has been used extensively to study focal traumatic brain injury. 
Although the impact variables can be well defined, little is known about the biomechanical trauma 
as delivered to different brain regions. The FE analysis based on high resolution T2- weighted MPJ 
images of rat brain was used in [Pena et al 20 05] to simulate displacement, mean stress, and shear 
stress of brain during impact. Young's Modulus E, to describe tissue elasticity, was assigned to each 
FE in three scenarios: in a constant fashion (E = 50 kPa), or according to the MRI intensity in a linear 
(E = [10, 100] kPa) and inverse-linear fashion (E = [100, 10] kPa). Simulated tissue displacement did 
not vary between the 3 scenarios, however mean stress and shear stress were largely different. The 
linear scenario showed the most likely distribution of stresses. 

A detailed FE model of the rat brain was developed in [Mao et al 2 006 for the prediction of 
intracranial responses due to different impact scenarios. The FE model was used to predict biome- 
chanical responses within the brain due to controlled cortical impacts (CCI). A total of six different 
series of CCI studies, four with unilateral craniotomy and two with bilateral craniotomy, were simu- 
lated and the results were systematically analyzed, including strain, strain rate and pressure within 
the rat brain. Simulation results indicated that intracranial strains best correlated with experimen- 
tally obtained injuries. An automating meshing method for patient-specific FE model was developed 
in [Guan et al 2006] . 3D geometries of two six-month-old infant heads were reconstructed from the 
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CT data. FE meshes including cranial bone of skull, brain, and suture were generated. Both static 
and dynamic analyzes were performed to verify the models. The study for blunt impact of infant head 
was performed by using these patient-specific models. 

A 3D FE analysis of human head was performed in Zo ng et al 2006| , to assess injury likelihood 
of the head subjected to impact loading. The structural intensity (SI) methodologjH was introduced 
in accordance with the prevailing practice in experimental biomechanics. The SI field inside the head 
model was computed for three frontal, rear and side impacts. The results for the three cases revealed 
that there existed power flow paths. The skull was, in general, a good energy flow channel. This study 
also revealed the high possibility of spinal cord injury due to wave motion inside the head. Recently, 
a 3D FE analysis was performed in Takahashi ct al 2 007] in respect to the frequency analysis of the 
pressure changes related to TBI. From the results of computer simulations and impact experiments, 
the authors found similar spectrums in some frequency bands, which indicated the occurrence of the 
brain injury. A vigorous shaking and an inflicted impact were compared in R,oth et al 200 7 , defined 
as the terminal portion of a vigorous shaking, using a FE model of a 6-month-old child head. Whereas 
the calculated values in terms of shearing stress and brain pressure remain different and corroborate 
the previous studies based on angular and linear velocity and acceleration, the calculated relative 
brain and skull motions that can be considered at the origin of a subdural haematoma show similar 
results for the two simulated events. A 2D FE model was developed in [Li et al 2 007] with objective 
to determine localized brain's strains in lateral impact using finite element modelling and evaluate the 
role of the falx. Motions and strains from the stress analysis matched well with experimental results 
from literature. A parametric study was conducted by introducing flexible falx in the finite element 
model. For the model with the rigid falx, high strains were concentrated in the corpus callosum, 
whereas for the model with the flexible falx, high strains extended into the cerebral vertex. 

On the other hand, it has been a common perception that rapid head rotation is a major cause of 
brain damage in automobile crashes and falls. A model for rotational acceleration about the center 
of mass of the rabbit head was presented in [Gutierrez et al 2001 , which allowed the study of brain 
injury without translational acceleration of the head. In the companion paper [Runncrstam et al 200 lj 
it was shown that rotational head acceleration caused extensive subarachnoid hemorrhage, focal tissue 
bleeding, reactive astrocytosis, and axonal damage. The initial response of the brain after rotational 
head injury involved brain edema after 24 h and an excitotoxic neuronal micro-environment in the first 
hour, which leaded to extensive delayed neuronal cell death by apoptosis necrosis in the cerebral cortex, 
hippocampus and cerebellum. Similarly, the study by | Zhang et al 2006] used the SIMon human 
FE head model and delineated the contributions of these accelerations using post mortem human 
subject (PMHS) lateral head impact experimental data. Results indicated that rotational acceleration 
contributed more than 90% of total strain, and translational acceleration produced minimal strain. 

The fidelity of cell culture TBI-simulations that yield tolerance and mechanistic information relies 
on both the cellular models and mechanical insult parameters. An electro-mechanical cell shearing 
device was designed by [LaPlaca et al 2005] in order to produce a controlled high strain rate injury 
(up to 0.50 strain, 30 s(-l) strain rate) that deforms 3D neural cultures (neurons or astrocytes in 
an extracellular matrix scaffold). Theoretical analysis revealed that these parameters generated a 
heterogeneous 3D strain field throughout the cultures that was dependent on initial cell orientation 
within the matrix, resulting in various combinations of normal and shear strain. 

Rigid-body modelling (RBM) was used in [Wolfson ct al 2005] to investigate the effect of neck 
stiffness on head motion and head-torso impacts as a possible mechanism of injury. Realistic shaking 
data obtained from an anthropometric test dummy (ATD) was used to simulate shaking. In each 
study injury levels for concussion were exceeded, though impact-type characteristics were required to 
do so in the neck stiffness study. Levels for the type of injury associated with the syndrome were not 
exceeded. 

The nonlinear mechanical behavior of porcine brain tissue in large shear deformations was de- 
termined in |Hrapko et al 2006] . An improved method for rotational shear experiments was used, 
producing an approximately homogeneous strain field and leading to an enhanced accuracy. The 
model was formulated in terms of a large strain viscoelastic framework and considers nonlinear vis- 
cous deformations in combination with non-linear elastic behavior. 

The relative motion of the brain with respect to the skull has been widely studied to investigate 

3 SI is a vector quantity indicating the direction and magnitude of power flow inside a dynamically loaded structure. 
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brain injury mechanisms under impacts, but the motion patterns are not yet thoroughly understood. 
The study of [Zou ct al 2007J analyzed brain motion patterns using the most recent and advanced 
experimental relative brain/skull motion data collected under low-severity impacts. With a minimum 
total pseudo-strain energy, the closed-form solutions for rigid body translation and rotation were 
obtained by matching measured neutral density target (NDT) positions with initial NDT positions. 
The brain motion was thus separated into rigid body displacement and deformation. The results 
showed that the brain had nearly pure rigid body displacement at low impact speed. As the impact 
became more severe, the increased brain motion primarily was due to deformation, while the rigid 
body displacement was limited in magnitude for both translation and rotation. Under low-severity 
impacts in the sagittal plane, the rigid body brain translation had a magnitude of 4-5 mm, and the 
whole brain rotation was on the order of +/-5 degrees. 

Biomechanical studies using postmortem human subjects (PMHS) in lateral impact have focused 
primarily on chest and pelvis injuries, mechanisms, tolerances, and comparison with side impact 



dummies. The objective of Yoganandan et al 2008 was to determine lateral impact-induced 3D 



temporal forces and moments at the head-neck junction and cranial linear and angular accelerations 
from sled tests using PMHS and compare with responses obtained from an anthropomorphic test 
device (dummy) designed for lateral impact. Results indicated that profiles of forces and moments at 
the head-neck junction and cranial accelerations were similar between the two models. However, peak 
forces and moments at the head-neck junction, as well as peak cranial linear and angular accelerations, 
were lower in the dummy than PMHS. Peak cranial angular accelerations were suggestive of mild TBI 
with potential for loss of consciousness. 

Deformation of the human brain was measured in vivo by [Sabet et al 2008] in tagged magnetic 
resonance images (MPJ) obtained dynamically during angular acceleration of the head, in order to 
provide quantitative experimental data to illuminate the mechanics of TBI. Mild angular acceleration 
was imparted to the skull of a human volunteer inside an MR scanner, using a custom MR-compatible 
device to constrain motion. Deformation of the brain was characterized quantitatively via Lagrangian 
strain. Consistent patterns of radial-circumferential shear strain occurred in the brain, similar to those 
observed in models of a viscoelastic gel cylinder subjected to angular acceleration. It has been noted, 
however, that strain fields in the brain are clearly mediated by the effects of heterogeneity, divisions 
between regions of the brain (such as the central fissure and central sulcus) and brain's tethering and 
suspension system, including the dura mater, falx cerebri, and tentorium membranes. 

The present paper proposes a new approach to brain injury dynamics, phrased as a coupled loading- 
rate hypothesis for TBI, stating that the main cause of TBI is an external Euclidean jolt, symbolically 
an Si? (3) —jolt, an impulsive loading striking the head in several degrees-of-freedom (both transla- 
tional and rotational) combined. This new concept is radically different from any of the standard 
FEM techniques proposed so far for brain injury mechanics (see the above literature review), as well 
as from any kind of Newton-Eulerian or Lagrangian/Hamiltonian injury dynamics, emphasizing 3 
new aspects of brain injury mechanics: (i) coupling of all 6 degrees-of-freedom; (ii) jolt dynamics 
rather than the force dynamics; and (iii) coupled dislocations/disclinations in the Cosserat multipo- 
lar viscoelastic continuum brain model. This new concept is a derivation of our previously defined 
concept of the covariant force law [Ivancevic and Ivancevic 2006bl llvancevic and Ivancevic 2006cl 
llvancevic and Ivancevic 20 07c . To support this hypothesis, we develop the coupled Newton-Euler 
dynamics of the brains's micro-motions within the cerebrospinal fluid (see Figure [T]), and from it de- 
rive the SE(3)— jolt dynamics, as well as its biophysical consequences in the form of brain's dislocations 
and disclinations. 



2 The 5E(3)-jolt: the cause of TBI 

In the language of modern dynamics |Ivancevic and Ivan cevic 20 06b| llvancevic and Ivancevic 20 06c, 
llvancevic and Ivancev ic 2007d, llvancevic and Ivancevic 2007ej . the microscopic motion of human brain 
within the skull is governed by the Euclidean SE(3)-group of 3D motions (see next subsection). Within 
brain's SE(3)-group we have both SE(3)-kinematics (consisting of SE(3)-velocity and its two time 
derivatives: SE(3)-acceleration and SE(3)-jerk) and SE(3)-dynamics (consisting of SE(3)-momentum 
and its two time derivatives: SE(3)-force and SE(3)-jolt), which is brain's kinematics x brain's mass- 
inertia distribution. 
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Informally, the external SE(3)-jol10 is a sharp and sudden change in the SE(3)-force acting on 
brain's mass-inertia distribution (given by brain's mass and inertia matrices). That is, a 'delta'-change 
in a 3D force-vector coupled to a 3D torque- vector, striking the head-shell with the brain immersed 
into the cerebrospinal fluid. In other words, the SE(3)-jolt is a sudden, sharp and discontinues shock 
in all 6 coupled dimensions of brain's continuous micro-motion within the cerebrospinal fluid (Figure 
[lj, namely within the three Cartesian (x, y, z)-translations and the three corresponding Euler angles 
around the Cartesian axes: roll, pitch and yaw. If the SE(3)-jolt produces a mild shock to the 
brain (e.g., strong head shake), it causes mild TBI, with temporary disabled associated sensory-motor 
and/or cognitive functions and affecting respiration and movement. If the SE(3)-jolt produces a hard 
shock (hitting the head with external mass), it causes severe TBI, with the total loss of gesture, speech 
and movement. 

The SE(3)-jolt is rigorously defined in terms of differential geometry [Ivancevic and Ivancevic 2006cl 
llvancevic and Ivancevic 2007ellTvancevic and Ivancevic 2007c| . Briefly, it is the absolute time-derivative 
of the covariant force 1-form (or, co- vector field). The fundamental law of biomechanics is the covari- 
ant force law [Ivancevic and Ivancevic 20 06b, Ivancevi c and Ivancevi c 2006c, Ivancevi c and Ivancevic 2007ej . 
which states: 

Force co-vector field = Mass distribution x Acceleration vector-field, 

which is formally written (using the Einstein summation convention, with indices labelling the three 
Cartesian translations and the three corresponding Euler angles): 

Fn = m^a" , (fi, v = 1, ...,6) 

where denotes the 6 covariant components of the external "pushing" SE(3)-force co-vector field, 
rripu represents the 6x6 covariant components of brain's inertia-metric tensor, while a v corresponds 
to the 6 contravariant components of brain's internal SE(3)-acceleration vector-field. 

Now, the covariant (absolute, Bianchi) time-derivative of the covariant SE(3)-force defines 
the corresponding external "striking" SE(3)-jolt co-vector field: 

§ t (F,) = m^|(o") = (a v + T^a x ) , (1) 

where ^(a") denotes the 6 contravariant components of brain's internal SE(3)-jerk vector-field and 
overdot (') denotes the time derivative. T^ x are the Christoffel's symbols of the Levi-Civita connection 
for the SE(3)-group, which are zero in case of pure Cartesian translations and nonzero in case of 
rotations as well as in the full-coupling of translations and rotations. 

In the following, we elaborate on the SE(3)-jolt concept (using vector and tensor methods) and 
its biophysical TBI consequences in the form of brain's dislocations and disclinations. 

2.1 S E(3) — group of brain's micro— motions within the CSF 

The brain and the CSF together exhibit periodic microscopic translational and rotational motion 
in a pulsatile fashion to and from the cranial cavity, in the frequency range of normal heart rate 
(with associated periodic squeezing of brain's ventricles) iMaier et al 1 9941. This micro-motion is 
mathematically defined by the Euclidean (gauge) SE(3) — group. Briefly, the SE(3) — group is defined 
as a semidirect (noncommutative) product of 3D rotations and 3D translations, 

SE(3) := SO(3) OR 3 . 

Its most important subgroups are the following (see Appendix for technical details): 

4 The mechanical SE(3)-jolt concept is based on the mathematical concept of higher— order tangency (rig- 
orously defined in terms of jet bundles of the head's configuration manifold) Ivancevic and Ivancevic 2006c 
Ivancevic and Ivancevic 2007c , as follows: When something hits the human head, or the head hits some external 
body, we have a collision. This is naturally described by the SE(3)-momentum, which is a nonlinear coupling of 3 linear 
Newtonian momenta with 3 angular Eulerian momenta. The tangent to the SE(3)— momentum, defined by the (absolute) 
time derivative, is the SE(3)— force. The second-order tangency is given by the SE(3)— jolt, which is the tangent to the 
SE(3)— force, also defined by the time derivative. 
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Subgroup 



Definition 



5*0(3), group of rotations 
in 3D (a spherical joint) 



Set of all proper orthogonal 
3 x 3 — rotational matrices 



SE(2), special Euclidean group 
in 2D (all planar motions) 



Set of all 3 x 3 — matrices: 
cos 8 sin 8 r x 
— sin 8 cos 8 r y 
1 



5*0(2), group of rotations in 2D 
subgroup of 5£?(2)-group 
(a revolute joint) 



Set of all proper orthogonal 
2 x 2 — rotational matrices 
included in SE(2) — group 



, group of translations in 3D 
(all spatial displacements) 



Euclidean 3D vector space 



In other words, the gauge SE(3) — group of Euclidean micro-motions of the brain immersed in 



the cerebrospinal fluid within the cranial cavity, contains matrices of the form 



R 





where 



b is brain's 3D micro-translation vector and R is brain's 3D rotation matrix, given by the product 
R = R v ■ R^ ■ Re of brain's three Eulerian micro-rotations, roll = R v , pitch = R^,, yaw = Rg, 
performed respectively about the a:— axis by an angle <£, about t he y— axis by an angle ip, and about 
the z— axis by an angle 8 |rvancevic""2 004, Park and Chung 2005] Hvancevic 2006] . 



R<n 



1 

cos ip — sin cp 
sin (p cos ip 



R* = 



cos ip sin ip 

1 
— sin ip cos ip 



Re = 



cos 9 — sin 9 
sin 9 cos 9 
1 



Therefore, brain's natural SE(3)— dynamics within the cerebrospinal fluid is given by the coupling 
of Newtonian (translational) and Eulerian (rotational) equations of micro-motion. 

2.2 Brain's natural SE(3)— dynamics 

To support our coupled loading-rate hypothesis, we formulate the coupled Newton-Euler dynamics of 
brain's micro- motions within the scull's SE(3)— group of motions. The forced Newton-Euler equations 
read in vector (boldface) form 



Newton 
Euler 



(2) 



p = Mv = F|pxu, 
TEl(i = T + 7rxw + pxv, 

where x denotes the vector cross productjfl 

M = My = diag{mi, m,2, m^} and I = = diag{I\, I2, -Z3}, (i, j = 1, 2, 3) 

are brain's (diagonal) mass and inertia matrices^ defining brain's mass-inertia distribution, with 
principal inertia moments given in Cartesian coordinates (x, y, z) by volume integrals 

I\= iff p{z 2 + y 2 )dxdydz, I2 = / / / p{x 2 + z 2 )dxdydz, I3 = / / / p(x 2 + y 2 )dxdydz, 



and 



dependent on brain's density p = p(x, y, z), 



5 Recall that the cross product u X v of two vectors u and v equals u X v = u-usinSn, where 6 is the angle between 
u and v, while n is a unit vector perpendicular to the plane of u and v such that u and v form a right-handed system. 

6 In reality, mass and inertia matrices (M.I) are not diagonal but rather full 3x3 positive-definite symmetric 
matrices with coupled mass- and inertia-products. Even more realistic, fully-coupled mass— inertial properties of a 
brain immersed in (incompressible, irrotational and inviscid) cerebrospinal fluid are defined by the single non-diagonal 
6x6 positive— definite symmetric mass-inertia matrix JVl se(3) i the so-called material metric tensor of the SE(3)— group, 
which has all nonzero mass— inertia coupling products. In other words, the 6x6 matrix A4gE(S) contains: (i) brain's 
own mass plus the added mass matrix associated with the fluid, (ii) brain's own inertia plus the added inertia matrix 
associated with the potential flow of the fluid, and (iii) all the coupling terms between linear and angular momenta. 
However, for simplicity, in this paper we shall consider only the simple case of two separate diagonal 3x3 matrices 
(M,I). 
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(where [ ]' denotes the vector transpose) are brain's linear and angular velocity vector^ (that is, 
column vectors), 

F = Fi = [FuFi, F 3 ] and T = T t = [T U T 2 , T 3 ] 

are gravitational and other external force and torque co-vectors (that is, row vectors) acting on the 
brain within the scull, 

p = p t . = Mv = \pi,P2,Pa] = [miVi,m 2 v 2 ,m 2 v 2 ] and 
7r = iTi = la; = [7ri,7T2,7r 3 ] = [h (Jith^h^a] 

are brain's linear and angular momentum co-vectors. 

In tensor form, the forced Newton-Euler equations ([2]) read 

Pi = MijiP = Fi + £ J ik PjUJ k , {i,j,k = 1,2,3) 
■hi = =^ + 4^^ + e 3 ikPj v k , 

where the permutation symbol e J ik is defined as 

f +1 if (i,j,k) is (1,2, 3), (3, 1,2) or (2,3,1), 
4 = < -1 if (ij,k) is (3, 2,1), (1,3, 2) or (2,1,3), 
I otherwise: i = j or j = k or k = i. 

In scalar form, the forced Newton-Euler equations ([2|) expand as 

f p 1 =Fi — m 3 v 3 uj 2 + m 2 v 2 uj 3 
Newton : I p 2 = F 2 + m 3 v 3 u>i — miv\uj 3 , (3) 



Pi 


= F 1 


- m 3 v 3 uj 2 + m 2 v 2 uj 3 




Po 


= F 2 


+ m 3 v 3 LUi - miviu 3 , 






= F 3 


- m 2 v 2 LUi + m 1 v 1 oj 2 




*1 


= Ti 


+ (m 2 - m 3 )v 2 v 3 + (I 2 


- I 3 )uj 2 uj 3 


7T 2 


= T 2 


+ (m 3 - mi)viv 3 + (I 3 




^3 


= T 3 


+ (mi - m 2 )viv 2 + (Ii 


- I 2 )loiuj 2 



Euler 

showing brain's individual mass and inertia couplings. 

Equations (HJ)-© can be derived from the translational + rotational kinetic energy of the braird 

E k = iv t Mv+-w*Iw, (4) 

or, in tensor form 

E = i-MyuV + -ItfwV. 

For this we use the Kirchhoff-Lagrangian equations (see, e.g., [Lamb 1932, Leonard 1997], or the 
original work of Kirchhoff in German) 

j t d v E k = 8 v E k xu + F, (5) 

id u E k = d u E k x u + d v E k x v + T, 
at 

where d v E k = d ul E k = in tensor form these equations read 

d 



-d vi E = e 3 ik (d vi E)w* + F i 

d_ 

db 



d^E = e° ik (d uj E)u J k +4 k (d vj E)v k + T i . 



7 In reality, u is a 3 X 3 attitude matrix (see Appendix). However, for simplicity, we will stick to the (mostly) 
symmetrical translation— rotation vector form. 

8 In a fully— coupled Newton-Euler brain dynamics, instead of equation we would have brain's kinetic energy 
defined by the inner product: 
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Using ([1]H([5]), brain's linear and angular momentum co- vectors are denned as 

p = d v E k , n = d u E k , 

or, in tensor form 

Pi = d v ,E, ir i =d ul .E, 
with their corresponding time derivatives, in vector form 

d d ^ v ■ d d a jt 

P= dt P= dt dvE > * =dt n =dt d " E > 

or, in tensor form 

d d a o ■ d d 

ft = dt Pi = dt d " iE > ni = dt ni = dt^ 

or, in scalar form 

P = [pi,P2)P3] = [mxiximiV^msbs], n = [tti, tt 2 , tt 3 ] = [I\Cj\, I 2 u>2, h^]- 

While brain's healthy SE( 3)— dynamics within the cerebrospinal fluid is given by the coupled 
Newton-Euler micro-dynamics, the TBI is actually caused by the sharp and discontinuous change 
in this natural SE(3) micro-dynamics, in the form of the SE(3) — jolt, causing brain's discontinuous 
deformations. 



2.3 Brain's traumatic dynamics: the SE(3) — jolt 

The SE(3)~ jolt, the actual cause of the TBI (in the form of the brain's plastic deformations), is 
defined as a coupled Newton+Euler jolt; in (co)vector form the SE(3) — jolt reada^l 



SE(3) - jolt 



Newton jolt :F = p- pxw-pxii , 

Eulcr jolt : T = 7r — -k xw-7rxu-pxv-pxv, 



where the linear and angular jolt co- vectors are 

F = Mv = [F U F 2 ,F 3 ], ± = lu=[f 1 ,f 2 ,t 3 ], 

where 

v=[vi,V2,V3]*, to = [0)1,0)2,0)3]', 

are linear and angular jerk vectors. 

In tensor form, the SE(3)— jolt readJ"1 

Fi = pi- e J lkP] uj k - s J ikPj uj k , k = 1, 2, 3) 

= - 4k*j" k - 4k^ k - 4kPi vk - 4 k Pjv k 

in which the linear and angular jolt covectors are defined as 

F = Fi = Mv = My i) J = [F u F a ,Fa], 
f = T i = K>=7 iJ <y I ' = [fi,f2,f3], 

where v = v i , and u) = uj 1 are linear and angular jerk vectors. 



9 Note that the derivative of the cross-product of two vectors follows the standard calculus product-rule: ^(uxv) = 
u x v + u x v. 

10 In this paragraph the overdots actually denote the absolute Bianchi (covariant) time-derivative JTJ, so that the jolts 
retain the proper covector character, which would be lost if ordinary time derivatives are used. However, for the sake 
of simplicity and wider readability, we stick to the same overdot notation. 
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In scalar form, the SE(3) — jolt expands as 

Pi - m 2 uu 3 v 2 + m 3 [uo 2 v 3 + v 3 cu 2 ) - m 2 v 2 6j 3 , 
p 2 + m\U)zV\ — m 3 uiv 3 - m 3 v 3 LUi + m\V\ib 3 , 
p 3 - miLU 2 vi + m 2 ujiv 2 - v 2 uJi - miviuj 2 , 

7fi - (m 2 - m 3 ) (v 3 v 2 + v 2 v 3 ) - (I 2 - I 3 ) {lo 3 lo 2 + U> 2 UJ 3 ) , 
ir 2 + (mi - m 3 ) (v 3 v-i + viii 3 ) + [h - I 3 ) [uj 3 6j 1 + u x u) Z ) , 
7r 3 - (mi - m 2 ) (w 2 wi + Vii) 2 ) - (Ii - I 2 ) {lo 2 lu 1 + u>iU) 2 ) ■ 

We remark here that the linear and angular momenta (p,7r), forces (F, T) and jolts (F,T) are 
co- vectors (row vectors), while the linear and angular velocities (v,u>), accelerations (v,u;) and jerks 
(v,£>) are vectors (column vectors). This bio-physically means that the 'jerk' vector should not be 
confused with the 'jolt' co-vector. For example, the 'jerk' means shaking the head's own mass-inertia 
matrices (mainly in the atlanto-occipital and atlanto-axial joints), while the 'jolt'means actually hit- 
ting the head with some external mass-inertia matrices included in the 'hitting' SE(3)-jolt, or hitting 
some external static/massive body with the head (e.g., the ground - gravitational effect, or the wall - 
inertial effect). Consequently, the mass- less 'jerk' vector represents a (translational+rotational) non- 
collision effect that can cause only weaker brain injuries, while the inertial 'jolt' co-vector represents 
a (translational+rotational) collision effect that can cause hard brain injuries. 

For example, while driving a car, the SE(3)-jerk of the head-neck system happens every time the 
driver brakes abruptly. On the other hand, the SE(3)-jolt means actual impact to the head. Similarly, 
the whiplash-jerk, caused by rear-end car collisions, is like a soft version of the high pitch-jolt caused 
by the boxing 'upper-cut'. Also, violently shaking the head left -right in the transverse plane is like a 
soft version of the high yaw-jolt caused by the boxing 'cross-cut'. 



Newton jolt 




Eulcr jolt 




2.4 Brain's dislocations and disclinations caused by the S E(3)— jolt 

Recall from introduction that for mild TBI, the best injury predictor is considered to be the prod- 
uct of brain's strain and strain rate, which is the standard isotropic viscoelastic continuum concept. 
To improve this standard concept, in this subsection, we consider human brain as a 3D anisotropic 
multipolar Cosserat viscoelastic continuum [Cosserat and Cosscrat 1898, Cosscr at and Cosserat 1 909, 
Eringcn 2002 1, exhibiting coupled-stress-strain elastic properties. This non-standard continuum model 
is suitable for analyzing plastic (irreversible) deformations and fracture mechanics Bilby and Eshclb y 1968| 
in multi-layered materials with microstructure (in which slips and bending of layers introduces ad- 
ditional degrees of freedom, non-existent in the standard continuum models; see jMindlin 19651 
Lak es 19'85] for physical characteristics and |Yang and Lakes 1981 Yang and Lakes 1982, Par k and Lakes 1 986] 
for biomechanical applications). 

The SE(3)— jolt (F,T) causes two types of brain's rapid discontinuous deformations: 

1. The Newton jolt F can cause micro-translational dislocations, or discontinuities in the Cosserat 
translations; 

2. The Euler jolt T can cause micro-rotational disclinations, or discontinuities in the Cosserat 
rotations. 

For general treatment on dislocations and disclinations related to asymmetric discontinuous defor- 
mations in multipolar materials, see, e.g., |Jian and Xiao-ling 1 995, Yan g et al 200 1| . 

To precisely define brain's dislocations and disclinations, caused by the Si? (3) —jolt (F,T), we 
first define the coordinate co- frame, i.e., the set of basis 1-forms {dx 1 }, given in local coordinates 
x l = (2; 1 , x 2 , a; 3 ) = (x,y,z), attached to brain's center-of-mass. Then, in the coordinate co-frame 
{dx 1 } we introduce the following set of brain's plastic-deformation-related SE(3)— based differential 
p— formal (see, e.g., [Ivancevic and Ivancevic 2006c| llvancevic and Ivancevic 2007e] ): 



11 Differential p— forms are totally skew-symmetric covariant tensors, defined using the exterior wedge-product and 
exterior derivative. The proper definition of exterior derivative d for a p— form (3 on a smooth manifold M, includes the 
Poincare lemma Ivancevic and Ivancevic 2006c, Ivancevic and Ivancevic 2007e : d(df3) = 0, and validates the general 



11 



the dislocation current 1-form, J = J i dx % ; 
the dislocation density 2-form, a = dx 1 A dx^\ 
the disclination current 2-form, S = \Sij dx 1 A dx^; and 
the disclination density 3-form, Q = -^Qij^dx 1 Adx^ A dx k , 
where A denotes the exterior wedge-product. According to Edelen Edelen 1980, Kadic and Edelen 1983 , 
these four SE(3)— based differential forms satisfy the following set of continuity equations: 

a = -dJ - S, (6) 
Q = -dS, (7) 
da = Q, (8) 
dQ = 0, (9) 

where d denotes the exterior derivative. 

In components, the simplest, fourth equation ©, representing the Bianchi identity, can be rewrit- 
ten as 

dQ = diQ[ ljk ] dx 1 A dx 1 A dx j A dx k = 0, 

where di = d/dx l , while Ouj...] denotes the skew-symmetric part of 0-y.... 
Similarly, the third equation (JSJ) in components reads 



Qijkdx 1 Adx j Adx k = dkOt[ij] dx k A dx 1 A dx\ 



1 
3! 

Qijk = -6dk<X[ij]- 
The second equation ([7]) in components reads 

T^Qijk dx 1 A dx j A dx k = -df.S[ij] dx k A dx 1 A dx\ 
Qijk — 6dkS[ifl. 
Finally, the first equation ([5]) in components reads 

-ciy dx 1 A dx-' = {djJi — -^Sij) dx 1 A dx-' , i 

In words, we have: 



The 2-form equation © defines the time derivative a =\oiij dx 1 Adx^ of the dislocation density 
a as the (negative) sum of the disclination current S and the curl of the dislocation current J. 



• The 3-form equation ([7]) states that the time derivative Q =-^Qijk dx 1 A dx^ A dx k of the discli- 
nation density Q is the (negative) divergence of the disclination current S. 



Stokes formula 

f P= f dp, 
JdM Jm 

where M is a p— dimensional manifold with a boundary and dM is its (p— 1)— dimensional boundary, while the integrals 
have appropriate dimensions. 

A p— form (3 is called closed if its exterior derivative is equal to zero, 

d/3 = 0. 

From this condition one can see that the closed form (the kernel of the exterior derivative operator d) is conserved 
quantity. Therefore, closed p— forms possess certain invariant properties, physically corresponding to the conservation 
laws. 

A p— form p that is an exterior derivative of some (p — 1)— form a, 

(3 = da, 

is called exact (the image of the exterior derivative operator d). By Poincare lemma, exact forms prove to be closed 
automatically, 

d/3 = d(da) = 0. 

This lemma is the foundation of the de Rham cohomology theory Ivancevic and Ivanccvic 2006c 
Ivanccvic and Ivancevic 2007c, Ivancevic and Ivanccvic 2007c . 
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• The 3-form equation JSJ) defines the disclination density Q as the divergence of the dislocation 
density a, that is, Q is the exact 3-form. 

• The Bianchi identity © follows from equation Q by Poincare lemma and states that the 
disclination density Q is conserved quantity, that is, Q is the closed 3-form. Also, every 4-form 
in 3D space is zero. 

From these equations, we can derive two important conclusions: 

1. Being the derivatives of the dislocations, brain's disclinations are higher-order tensors, and thus 
more complex quantities, which means that they present a higher risk for the severe TBI than 
dislocations — a fact which is supported by the literature (see review of existing TBI-models 
given in Introduction). 

2. Brain's dislocations and disclinations are mutually coupled by the underlaying SE(3)— group, 
which means that we cannot separately analyze translational and rotational TBIs — a fact which 
is not supported by the literature. 

3 Conclusion 

Based on the previously developed covariant force law, in this paper we have formulated a new 
coupled loading-rate hypothesis for the TBI, which states that the main cause of traumatic brain 
injury is an external SE(3) — jolt, an impulsive loading striking the head in several degrees-of-freedom, 
both rotational and translational, combined^] To demonstrate this, we have developed the vector 
Newton-Euler mechanics on the Euclidean SE(3) — group of brain's micro- motions within the cere- 
brospinal fluid. In this way, we have precisely defined the concept of the SE(3) — jolt, which is a cause 
of brain's rapid discontinuous deformations: (i) translational dislocations, and (ii) rotational discli- 
nations. Based on the presented model, we argue that: (1) rapid discontinuous rotations present a 
higher risk for the severe TBI than rapid discontinuous translations, and (2) that we cannot separately 
analyze rapid brain's rotations from translations, as they are in reality coupled. 

4 Appendix: The S E(3)— group 

Special Euclidean group SE(3) :— 50(3) t> R 3 , (the semidirect product of the group of rotations with 
the corresponding group of translations), is the Lie group consisting of isometrics of the Euclidean 3D 
space R 3 . 

An element of SE{3) is a pair (A, a) where A e 50(3) and ael 3 . The action of SE(3) on R 3 is 
the rotation A followed by translation by the vector a and has the expression 

{A, a) ■ x = Ax + a. 

The Lie algebra of the Euclidean group SE(3) is se(3) =l 3 xK 3 with the Lie bracket 

[(£, u), (rj 1 v)} = {£xri^xv-rix u). (10) 

12 One practical application of the proposed model is in design of helmets. Briefly, a 'hard' helmet saves the skull 
but not the brain; alternatively, a 'soft' helmet protects the brain from the collision jolt but does not protect the skull. 
A good helmet is both 'hard' and 'soft'. In other words, if a human head covered with a solid helmet collides with a 
massive external body, the skull will be protected by the helmet — but the brain will still be shocked by the SE(3)-jolt, 
and a TBI will be caused. With or without the 'hard' helmet, brain's inertia tensor will be moved and rotated by the 
external SE(3)-jolt, and this will cause a brain injury, proportional to the jolt— collision with the head. Therefore, while 
protecting the skull is a necessary condition for protecting the brain, it is not enough. Brain's inertia tensor needs 
another kind of protection from the external collision-jolt. Contrastingly, if a human head covered with a 'soft' helmet 
collides with a massive external body, the helmet will dissipate the energy from the collision jolt, but will not necessarily 
protect the skull. As a result, a proper helmet would have to have both a hard external shell (to protect the skull) and 
a soft internal part (that will dissipate the energy from the collision jolt by its own destruction, in the same way as 
a car saves its passengers from the collision jolt by its own destruction). Note that, hypothetically speaking, an ideal 
shock— absorber is not a classical spring-damper system (with the distance-dependent spring and velocity-dependent 
damper), but rather a constant-resistance damper. 
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Using homogeneous coordinates, we can represent SE(3) as follows, 

SE(3) = ^ GGi(4,R):i?e50(3),peR 3 |, 

with the action on R 3 given by the usual matrix-vector product when we identify R 3 with the section 
R 3 x {1} c R 4 . In particular, given 

9 = ( J \ ) G S25(3), 

and g e R 3 , we have 

g ■ q = Rq + p, 



or as a matrix- vector product, 



The Lie algebra of SE(3), denoted se(3), is given by 

se(3) = {(jj ) G M 4 (R) : w £ so(3), v 6 IE 
where the attitude (or, angular velocity) matrix w : R 3 — > so (3) is given by 

























The exponential map, exp : se(3) — > SE{3), is given by 

exp 

where 












, 





exp(w) Av 
1 



1 — cos H^ll [Hjj-smjHI 2 

A — 1 i ^ W i 5 W , 

Ikll 2 P\f 

and exp(u>) is given by the Rodriguez' formula, 

, , r sinllwll 1 — cosllwll 2 

exp w = I + n I, ^ H ^7— . 

!MI |M| 
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